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Abstract 

We show that the multipulse application can suppress the decoherence of quantum entanglement 
with focusing on the concurrence, the degree of entanglement. By evaluating the time evolution of 
concurrence with a linearly interacting spin-boson model under pulse application, we find that the 
effectiveness of the multipulse control depends on the non-Markovian nature of the reservoir. 
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The application of quantum principle to the field of information processing has opened 
a new perspective for secure communication and high-speed computation. However, when 
we execute the schemes of quantum information processing, we have to overcome several 
obstacles which stem from the frailty of quantum substance. Above all, the vulnerability of 
quantum entanglement to a noisy reservoir constitutes a serious obstacle to the realization 
of quantum teleportationQ and quantum computation^]. 

Since the methods to combat against the decoherence of the entanglement have been 
firstly proposed for the schemes of quantum teleportationQ and quantum cryptography ^, 
many studies have been made under the condition/assumption that the entangled qubits are 
separated in space, and that we can only execute Local Operation and Classical Communi- 
cation(LOCC) after the qubit has been contaminated by a noisy reservoir. These studies, 
called as entanglement concentration, purification, or distillation 0, 0, Q, have tried to 
extract a smaller number of purely entangled pairs from two or more contaminated en- 
tangled pairs bv applying LOCC. While previous experiments have realized the scheme of 
distillation 0, Q or purification it is necessary to use high accuracy measurement of a 
qubit (photon) for classical communication and an almost infinite number of contaminated 
pairs of qubits to obtain a perfectly pure entangled pair; a process which seems practically 
unfeasible. 

In the field of quantum computation, many methods have been proposed to protect a 
single qubit against the effect of a noisy reservoir, which is called as decoherence. These are 
roughly categorized into three types: (1) methods to obtain one stable qubit with ancillary 



qubits 



12l . Il3| |. (2) methods to suppress the decoherence with quantum feed )ack[l^, and 



(3) application of pulse train (bang-bang or dynamical decoupling method) |15|, |16|. Above 
all, the bang-bang or dynamical decoupling method has attracted much attention for its 



feasibility to various kinds of svstems|l7||. However, the physical background of the bang- 
bang method was not fully discussed in these studies. In our recent work[3,Q, we discussed 
the effectiveness of the bang-bang method on the decoherence of a single qubit and showed 
that the key is in the partial reversibility of qubit, which arises from a non-Mar kovian nature 
of decoherence process. 

When a qubit interacts with a reservoir that has a finite correlation time, the phase relax- 
ation is described by an integral-differential equation for the off-diagonal element of density 
operator (or induced moment ) P{t) of qubit as ^P(t) = — dt' M(t')P(t — t') which is 
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obtained by the projection operator method j^]. The memory kernel M{t') is described by 
correlation functions of reservoir and often has significant values only when time t' is shorter 
than the correlation time of reservoir t^. This indicates that the time derivative of P{t) de- 
pends on its past history, which we call non-Markovian nature. When time goes much longer 
than Tc, the integral-differential equation is approximated as ^P{t) = —P{t)/T2 to obtain 
the exponential dephasing P{t) = e~*^^'^P{0) where T2 is defined as = dt' M{t'). 
The exponential decay, which is obtained by the above-mentioned Markovian approxima- 
tion, indicates that the time evolution is irreversible. Since the correlation time of reservoir 
Tc is usually much shorter than T2, the exponential dephasing is observed in various ma- 
terials. However, it should be noted that in the time scale which is comparable or shorter 
than Tc, P{t) shows the non-exponential decay, which cannot be described by the constant 
dephasing time T2. In other words, during such time scale, the time evolution of P{t) is 
partially reversible. The possibility to retrieve the reversibility from an induced moment 
with partial reversibility was firstly inv estig ated in the scheme of optical transient four- 
wave mixing (optical parametric effect) 21, 22, 2^ where application of a vr-pulse on 



an induced moment causes time-reversal action, and therefore P{t) is recovered depending 
on the degree of time-reversibility of the system. This was verified by the experiments for 
sodium resorufin in dimethylsulfoxide 125^ , iron-free myoglobin 26 1 , dye molecules j23], and 
CdSe quantum dots [2^. 

The bang-bang method using the vr-pulse train can be understood as a natural extension 
of the two pulse optical transients mentioned above, as was pointed out in [l^. That is, 
the repeatedly applied tt pulses with pulse interval much shorter than can prevent the 
decoherence of qubit because of the non-Markovian nature of the system. In the limit of 
zero pulse intervals in bang-bang method, the control of phase relaxation is possible even 
for the exponential dephasing. However, it should be noted that we should take the zero 
limit of pulse interval Tg before we take the zero limit of Tc; equivalently, the condition of 
Ts/tc should be satisfied when we take the limit ^ Therefore, Tc should be 
comparable or larger than the pulse interval in any case, and the non-Markovian property 
of decoherence of qubit is essential for the effectiveness of bang-bang control. 

In this paper, we show that the multipulse application can control entanglement for a pair 
of qubits. In order to analyze the controllability of entanglement by multipulse application, 
it is plausible to evaluate a quantity called concurrence [2^ that plays a center role to 
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describe the degree of entanglement of a pair of qubits. For a density matrix p of a pair of 
qubits under consideration, the concurrence C{p) is defined by 

C(p) =max(0, 2 A^ax-Tri?). (1) 
Here Amax is the maximum eigenvalue of the operator R which is defined by 

R = sf^Tp. (2) 
for the density matrix of the pair of qubits p and p, 

P = [Ox^y ® cr2,j/)P*(o-l,j/ ® Cr2,2/), (3) 

where the On^y indicates ?/-component of the Pauli matrix of n-th qubit and p* means the 
complex conjugate of density matrix p. Since the product of o"„ ^ and complex-conjugated 
operator is associated with time-reversal operator jSO'l, the concurrence is based on the "de- 

nn 

gree of equality" |22l, between a density matrix under consideration and a density matrix 
obtained by time-reversal operation. By this definition, we can consider that the concur- 
rence describes the degree of reversibility. This is exactly the reason why the concurrence 
is an appropriate quantity when we consider the vr pulse application to a pair of qubits 
that interact with a reservoir of non-Markovian nature, since the pulse application causes 
time-reversal in the evolution of the qubits. 

Taking a simple case where a pair of qubits linearly interact with a common reservoir, 
let us show the effectiveness of multipulse application to suppress the disentanglement. The 
Hamiltonian of this system is 

'Hr = 7^0 + 'HsB = CHs + Hb) + 'HsB , (4) 
with Hs = h uJoSn,z, Hb = Y.k h^kh\hk, and 

2 

UsB = h'^Sn,z'^hkUJk{hk + h\). (5) 

n=l k 

Here Sn,z is the z-component of the n-th qubit (n = 1,2); hk{h\) indicates the annihila- 
tion (creation) operator of the fc-th boson which composes the reservoir; hk is the coupling 
strength between the qubit; the fc-th boson of the reservoir. It is noted that the reser- 
voir is common for two qubits, and the model is applicable to a pair of quantum dots 
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in semiconductors [32] • To suppress the disentanglement due to the interaction with the 
reservoir, we apply sufficiently short 7r-pulse train on both qubits simultaneously. The 
Hamiltonian under pulse application is written as 

N 

-Hspit) = Us + Y^np.it), (6) 

i=o 

1 ^ 

Hp,j{t) = --E,{t) ■ Pi 5^(5„,+e-*-°* + ^„,_e*""*) (7) 

n=l 

where Ej{t) denotes the field amplitude of j-th applied pulse. We assume that two qubits 
have the same transition moment jl. When we consider a \ spin (a two-level system) as a 
qubit, Eq.® means the magnetic interaction (the electric interaction for optical transition), 
respectively . 

When we apply vr-pulses with pulse interval and pulse duration At, the density 
operator for the total system p{t) is written as 

N 

x{Y[e-'^^'r:-At'it'LsPit')^-iLnx{rs-At)^p^Q^^ (8) 

i=i 

where p(0) indicates the density operator for the total system at an initial time t = 0, 
and Lji and Lspit) are Liouville operators defined by L^X = ^[Hr,X] and Lsp(t)X = 
^[Hsp(t), X] for an arbitrary operator X. The exponential of Liouville operators are related 
to the exponentials of the corresponding Hamiltonians as e~*^^*X = e~fi^^*Xe'i^^* and 
^iI^::_^,d,'Ls,(^)^ = e-^^--A.rf*'{>^s+Wp.(*'))xei/^^-A.'^*'(^^+w-..(*')) ^here we assumed 

that the pulses are well separated. 

In order to evaluate the time evolution of the density operator, it is convenient to use a 
relation for the time evolution operator as e~^^^* = e~^^°*T+ exp[— | jl'HsB{i')dt'] where 
T+ denotes the time ordering symbol from right to left, and HsB{t) = e'f'^"^HsBe~^''^°*- 
When we consider the decoherence caused by Eq.(jSI), the time ordered exponential is written 
as 

i /■* ~ 

T+exp[-- nsB{t')dt'] = Diag[n+(t), 1, l,M_(t)], (9) 

where Diag[. . .] means a diagonal matrix with elements in the bracket [■■■], and u±{t) = 
T+exp[T| jQB{t')dt'] with B{t) = hY^f^hkUkipke-'"^"^ + fo^e^"^**). Assuming that we apply 
square pulses with height Ej for the j-th pulse, we obtain the time evolution operator 



under pulse application as e r //T^-At^^^'C^s+Hp.jCi')) ^ ^-j^nsAt^-j^np^jAt^ where Hpj means 
the pulse Hamiltonian in the interaction picture with Ti-pj = e^^^^'Hpj{t)e~i''^^* = —\Ej ■ 

Let us consider the following initial conditions for the qubits and the boson reservoir: 
(1) the two qubits maximally entangled at an initial time as = -^(|1)|1) + |0)|0)) = 
^(|11) + |00)) where |0) and |1) indicate the two states of a qubit, and (2) the boson 
reservoir is in the vacuum state. Then we have the density operator at the initial time as 

f PR pr\ 


\Pb O prJ 

where pp = |0)(0| is the density operator for the boson reservoir. 

Using Eq.(P3) ~ Eq.(|10p. we obtain the time evolution of concurrence under vr-pulses 
with pulse interval Tg in the limit of infinitely short pulse duration At — > as follows: 



p(o) = \ij){ij\^PR 



(10) 



C{t) = lexp[-2Y,\Mt)\' 



(11) 



where 



N 
m=l 



(12) 



Equation (lllj) is similar to the formula that we have obtained for the decoherence of 
a single qubit under linear interaction with the boson reservoir (see Eq.(29) in 2). The 
exponent in Eq.((7j) is twice larger than that in Eq.(29) in which is the signal intensity 
expressed by the off-diagonal element of the reduced density matrix as I(t) = \ (e|Trijp(t) |e) p 
where \e){\g)) corresponds to |0)(|1)), respectively. This means that the entanglement of 
two qubits decays four times faster than the single qubit, since the concurrence is given by 
the off-diagonal element of density matrix of two qubits as |(00|Tr/i;p(t)|ll)p. The same 
result has been pointed out by Yu and Eberly,.33]- 

To evaluate Eq. (jll|) . we need an explicit form of the coupling function 



h{u;) = J2\hkM 



(13) 



Bang-bang control is effective for any type of coupling function when the pulse train is 
applied with infinitely small pulse interval, which is practically difficult to execute. As we 
discussed for the decoherence of a single qubit jl^ . we can release such a difficulty by paying 
attention to the dynamical behavior of reservoir: by synchronizing the pulse-train applica- 
tion with the characteristic oscillation period of the reservoir, we can efficiently suppress the 
decoherence of the qubit, a process which we call "synchronized pulse control" (SPC)^^. 
The effectiveness of pulse control with finite pulse interval, especially for SPC, strongly 
depends on the type of coupling function. Fixing the width of the coupling function at a 
relatively small value, we compared the effectiveness for Lorentzian and non-Lorentzian cou- 
pling function. We found that SPC is more effective for non-Lorentzian coupling function 
than Lorentzian coupling function. We can explain the reason with a picture of two-step 
structured reservoir which is transformed from the original spin-boson model. When the 
coupling function is Lorentzian, we can reconstruct the reservoir which consists of infinite 
number of bosons to obtain a single boson, called an interaction modej^J] or quasi mode[3|, 
which interacts with another reservoir of infinite number of bosons with a flat (white) cou- 
pling function. In this case, the interaction mode shows the Markovian time evolution. 
When a single qubit interacts with such a boson with Markovian nature, the qubit loses 
partial reversibility quickly, which makes the pulse control less effective. Since Eq. fjllj) is 
similar to the time evolution of the decoherence of a single qubit under pulse application, 
we can say that the pulse control on entanglement is also less effective for Lorentzian cou- 
phng function. To discuss the effectiveness of pulse control with finite pulse interval, let us 
consider the non-Lorentzian coupling function. 

As a typical example of non-Lorentzian coupling function, we assume the coupling func- 
tion to be a Gaussian distribution with the mean frequency Up and the variance 7^, 



where s is the average number of bosons interacting with a qubit, which means the strength 
of coupling between a qubit and reservoir. 

In Fig.l, we show the time dependence of the concurrence wherein a scaled time 

variable t = tUp is used, and the parameters are set as 7^ = ^pjojp = 0.1 and s = 5. Without 
pulse application, one can see that the concurrence shows a damped oscillation with the 
center frequency ojp of the couphng function (the mean period is In in the scaled time) as 




(14) 
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shown in Fig. 1(a). In the figure, we can see that the partial reversibihty of concurrence 
remains during the mean period of dynamical motion of boson reservoir. When one applies 
a vr-pulse train with a relatively short interval = 7r/5, one can see that the damping of 
the concurrence is reduced and approaches to a constant value(Fig.l(b)). This result clearly 
indicates that the degradation of concurrence is effectively suppressed by the application of 
vr-pulse train. With increasing the pulse interval, the concurrence decreases. At the pulse 
interval fs = vr, the concurrence takes even smaller values than those in the case of no pulse 
control as shown in Fig. 1(c). However, when the pulse application is synchronized with the 
oscillation period of the center frequency of the coupling function by setting = 271, one 
obtains a periodically recovering concurrence. The recovery means that the time-reversal 
action by vr-pulse application efficiently recover the coherence of the quantum entanglement 
by paying attention to the dynamical motion of the boson reservoir which causes the memory 
(non-Markovian) effect on the time evolution of concurrence. 
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FIG. 1: Time evolution of C{t) for 7^ = ^p/ujp = 0.1, s = 5. (a) shows the case without pulse 
application, (b) for pulse interval = vr/S, (c) for = vr, (d) for fg = 2tt. 

As a measure of entanglement, we use the concurrence in this paper. When the interaction 
between the qubits and the reservoir is written as Eq.(j21), we can relate other measures with 
concurrence. Firstly, we consider the entropy of the pair of qubits S{t) = — J2t=i log4 
which has been proposed by Thew and Munro [s^ . They categorized the function of entan- 
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glement distillation, purification, and concentration by the dependence of the entanglement 
of formation(EOF) on the entropy S{t). They defined the concentration as the procedure 
to decrease the entropy and increase EOF, purification as that to decrease entropy without 
changing EOF, and distillation as that to increase EOF without changing entropy. For the 
present model where the system-reservoir interaction is adiabatic, we obtain a simple rela- 
tion between the time evolution of entropy and the concurrence as S(t) = — J2t=i -^j log4 -^i 
where Ai = _ ^ i-c'(f) -^2^ ^3 = = Q. Since this relation holds under the pulse 

excitation, the concurrence and the entropy S(t) are mutually related during the course of 
entanglement recovery by multipulse application, as Fig. 2 shows. The entropy S(t) ranges 
between and 0.5 in this case. This is because the density matrix of the most phase de- 
cohered state has only two elements as (00|Tr/jp(t)|00) = (ll|Tr/jp(t)|ll) = i . We also 
found the relation between the purity of qubits P{t) = Tr{p{ty') and the concurrence as 
P{t) = ^^^^j^, which shows how the multipulse application can recover the purity of the 
qubits as the concurrence decreases. 
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FIG. 2: The dependence of C{t) on S{t). 



We also consider the case in which each qubit individually interacts with its own boson 
reservoir: 

2 

nsB = hJ2 Sn,. Yl hu,.u^Sht} + 6if ). (15) 

n=l kn 

Here 6^."^ (^i"^^) indicates the creation (annihilation) operator of boson reservoir with which 
n-th qubit interacts. In Eq.()15p. UO}.^ is the frequency of the A;-th boson of the n-th reservoir. 

indicates the coupling strength between the n-th qubit and the A;-th boson. At an initial 
time, we assume the qubits to be in the maximally entangled state as before. In this case, 
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we obtain the expression for the time evolution of concurrence as, 

C^W = ^flexp[-l5^|a,„(t)|2], (16) 

n=l kn 

where afc„(^) is given by replacing suffixes k in Eq. (jl3p with k^- When we set the same 
coupling function for two reservoirs, InC(t) in non-common reservoir case is a half of that 
in the common reservoir case in Eq. pi|) . This means that the qualitatively the same effect 
of a TT-pulse train is observed for non-common reservoir case. Consequently, InC(t) becomes 
a half of that in the case of common reservoir shown in Eq.(|7j), if we set the same coupling 
function for two reservoirs. This is the only difference between the two cases of common 
and non-common reservoirs. Qualitatively the same effect of a vr-pulse train is observed. 

In this paper, we have analyzed the case where the two qubits do not interact each other 
directly, although the interaction between qubits is necessary to generate entanglement. 
This is because, when we assume the initial state of the pair of qubits to be = + 
1 00)), the Heisenberg interaction as S" ■ does not affect the time evolution of concurrence. 
Then, the interaction with common or non-common reservoir is sufficient to discuss the 
disentanglement in this case. 

In the evaluation of concurrence, we have used the Gaussian coupling function as an 
example. Since the formula of concurrence is found to be similar to the formula for a 
single qubit 1^] , we can say that the effectiveness of pulse control on the concurrence is also 
more effective for semi-elliptic coupling function than the Gaussian one, especially for the 
synchronized pulse control(SPC). This is because the semi-elliptic function does not have a 
tail, and is very different from the Lorentzian shape. However, when the pulse interval is 
infinitely small, the bang-bang control is effective for any type of coupling function. 

In summary, we have shown that the multipulse application can suppress the decoherence 
of quantum entanglement with focusing on the concurrence in the non-common reservoir case 
as well as in the common reservoir case. We found that the effectiveness of multipulse control 
depends on the non-Markovian nature of the reservoir. 
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